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Abstract. In the context of bubble universes produced by a first-order phase 
transition with large nucleation rates compared to the inverse dynamical time scale 
of the parent bubble, we extend the usual analysis to non-vacuum backgrounds. In 
particular, we provide semi-analytic and numerical results for the modified nucleation 
rate in FLRW backgrounds, as well as a parameter study of bubble walls propagating 
into inhomogeneous (LTB) or FLRW spacetimes, both in the thin-wall approximation. 
We show that in our model, matter in the background often prevents bubbles from 
successful expansion and forces them to collapse. For cases where they do expand, we 
give arguments why the effects on the interior spacetime are small for a wide range of 
reasonable parameters and discuss the limitations of the employed approximations. 
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1. Introduction 

In the recent literature a lot of interest has been dedicated to the question of how 
inflationary models [1, 2] can be embedded into a general theory [3]. In this respect 
the emergence of the string theory landscape [4] has opened various possibilities and 
a completely new perspective. Due to flux compactification in string theory there can 
be a huge number of distinct vacua [5] in scalar field space without a unique physical 
selection mechanism available at present. 

Classically, a field that is stuck in one vacuum would be trapped forever and could 
never move through the landscape. Quantum mechanically, the vacua are rendered 
metastable by the possibility for any field to tunnel to another (metastable) vacuum and 
thereby probe the landscape. The rate of tunneling naturally depends on the energy 
difference between the vacua and on the height and width of the potential barrier. The 
transition by pure tunneling can be described by the Coleman-De Luccia instanton [6] 
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Figure 1. We investigate the effects on the tunneling and the evolution of bubbles 
of new vacuum when the precursor state is non-vacuum. The simplistic sketch shows 
such a tunneling event that can be imagined in the context of rapid tunneling in the 
landscape. We are interested in the case when tunneling is rapid enough such that 
the background had no time to relax to vacuum and is, for instance, undergoing a 
cosmological phase transition. This complication can be relevant in the context of 
chain inflation [10]-[18] or other models that invoke rapid tunneling, like e.g. DBI 
or resonance tunneling [19], [20]. For some particular examples of non- vacuum 
backgrounds (like e.g. radiation dominated FLRW or inhomogeneous LTB), we analyse 
the effects on the semiclassical tunneling rate as well as on the subsequent general 
relativistic evolution of the formed new vacuum bubbles. 



whereas a transition over the barrier due to thermal fluctuations can be described by 
the Hawking-Moss instanton [7]. 

The Coleman-De Luccia process gives rise to the nucleation of spherically symmetric 
regions (bubbles) in space which are filled with new vacuum and expand into the 
old vacuum - a first-order phase transition. So, in the landscape multiple tunneling 
from multiple metastable vacua can occur, leading to different patches of spacetime 
undergoing coeval inflation in a variety of vacua. This process leads to a very 
complicated (fractal) large-scale structure of the universe that has been termed the 
attractor of eternal inflation. It is an attractor in the sense that statistical properties of 
the large-scale universe asymptotically do not depend on initial conditions [8]. 

The probability per unit four-volume for a tunneling to occur may be written as 

r = Aexp(-2Im5) , (1) 

where A is a determinantal factor and S is the action of the instanton mediating the 
tunneling process. The prefactor A is notoriously hard to obtain and is usually assumed 
to be of order unity for most practical purposes. The instanton action S can be obtained 
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by solving the field equations with appropriate boundary conditions. In many cases the 
layer separating the two phases can be thought of as a domain wall with a thickness 
that is small compared to the size of the bubble ('thin-wall approximation'). With this 
assumption a calculation in Minkowski spacetime [9] yields ImS = 277r^cr^/4e'^ , with e 
the latent heat and a the energy density in the surface layer. It is often assumed that 
cr ^ e (in Planck units) and thus tunneling is suppressed by a huge exponential factor, 
which implies extremely long lifetimes of the metastable vacua. However, any given 
point will eventually enter the decay chain, no matter what the odds are. 

In this work we will study bubbles which nucleate from a non- vacuum initial state, 
that is inside pockets which do not obey the de Sitter symmetries because they are 
not vacuum dominated at the time when tunneling takes place, see Fig. 1. We look at 
tunneling from a precursor state that is, say, in a radiation or a matter dominated phase 
(but not yet vacuum dominated), and which produces an inflationary bubble of de Sitter 
vacuum. One may object that a large class [21, 22, 23, 24] of inflating spacetimes, if 
they include a positive cosmo logical constant, will always be vacuum dominated in 
the asymptotic future (the cosmic no-hair conjecture [7, 25, 26]). Together with the 
aforementioned long lifetimes of metastable vacua in the landscape this seems to render 
the question of bubble nucleation from anything other than the vacuum academic. 
However, there are recent proposals of inflation where tunneling between minima on 
the landscape occurs rapidly. This is the case, e.g., in the model of chain inflation^ 
[10, 11, 12, 13, 14, 15, 16, 17, 18] or for non-standard tunneling on the landscape via 
DBI or resonance tunneling [19], [20]. In case of rapid tunneling - in chain inflation 
every transition only yields a fraction of an e-fold - the efficiency of the cosmic no-hair 
mechanism can be questioned. Therefore, we argue that the question of the consequences 
of rapid inflation through a pocket that has possibly not yet relaxed to pure vacuum 
deserves some attention. We will attack this problem in a twofold way. On the one 
hand, we will look at the tunneling itself and on the other hand, we will study the 
classical evolution of a vacuum bubble in a non- vacuum environment. 

Concerning the former, one would for instance expect that the instanton action 
picks up some geometrical corrections due to the non-trivial background. Any such 
modiflcations are potentially interesting, since the tunneling rates are exponentially 
sensitive to them. We will show how one can assess the relevance of those effects, based 
on a comparison of characteristic time scales. Geometrical corrections become important 
if some dynamical time scale of the background is comparable to or smaller than the 

I In chain inflation, tunneling occurs stepwise through a sequence of many minima. Chain inflation 
resurrects the old inflationary idea of a first-order phase transition but is able to solve the problem of 
graceful exit. Originally, this was accomplished because the fields were assumed to be coupled. The 
coupling is responsible for rapid tunneling: once a first tunneling has occurred the field increases the 
decay probability of its neighbor (s) further up in the landscape due to coupling. A chain reaction of 
rapid tunneling starts that eventually ends with a (slow) transition to true vacuum. However, chain 
infiation does not seem to allow for eternal inflation because it would not produce the right primordial 
density fluctuations [12]. In [13] a concrete realization of chain inflation on the string landscape that 
is driven by four form fluxes has been proposed. 
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characteristic time scale of the tunnehng process, which is given by the bubble's proper 
nucleation radius. This may happen either when the nucleation radius is very large - but 
then tunneling rates are so small that the background geometry will relax to a de Sitter 
stage long before transition occurs - or when the background dynamics is characterized 
by some very short time scales. 

The subsequent classical evolution of bubbles in such an environment will also 
be modified with respect to the vacuum case. We will study the propagation of 
vacuum bubbles in presence of i) homogeneously distributed matter, ii) matter with an 
inhomogeneous radial profile and iii) a fluid undergoing a second-order phase transition. 
These simple toy scenarios shall give a taste of the phenomenology which may result 
from non- vacuum initial states. 

Besides for chain inflation, the results of this analysis can be important also for other 
scenarios. For instance, the influence of background inhomogeneity can be relevant in 
the context of landscape sampling by tunneling. While through resonant processes the 
tunneling rate on the landscape can be enhanced [20], inhomogeneous initial states may 
be of high importance for this sort of tunneling [27, 28]. The approach in our work 
is partly inspired by [29], where it was studied how an inhomogeneous and spherically 
symmetric background influences the classical evolution of an inflationary bubble of new 
vacuum immersed into it. In the limits of the used model it was claimed in [29] that 
ambient inhomogeneities do not alter the evolution of bubbles significantly as long as 
the weak energy condition is respected. On the other hand, this topic also touches the 
interesting issue of inhomogeneous initial conditions versus the onset of inflation, see 
e.g. [30, 31, 32, 33, 34, 35]. 

The paper is organized as follows. In section 2 we present results that show the 
influence of dynamical backgrounds on the nucleation rate of bubbles of new vacuum. 
To this end, we apply an extension of the usual semiclassical approach to cosmologically 
interesting time-dependent FLRW backgrounds such as power law inflation or radiation 
dominated universes. The nucleation rates are obtained in the thin-wall approximation 
by using a complex time path formalism. In section 3, we analyze the subsequent 
classical trajectory of bubbles. In addition to FLRW backgrounds§, we look at the 
evolution in an exact spherically symmetric and inhomogeneous spacetime containing 
also matter. The interior of the bubble is assumed to be de Sitter spacetime, as a 
first approximation to the inflationary phase of the patch of the universe that we 
are observing. The bubble evolution follows from the Israel junction method which 
is employed to join these spacetimes together. We ask whether signatures of the 
background can potentially be seen by an observer inside the bubble. In section 4 
we summarize our results and give conclusions as well as some remarks on the limits of 
the used methods and an outlook. We use units c = h = G = 1, and sign conventions 
for geometrical quantities in accordance with [36]. 

§ We thank Ben Freivogel for this suggestion. 
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2. Bubble nucleation in time-dependent settings 

Quantum nucleation of bubbles is a very intricate problem, especially when effects of 
gravity have to be taken into account. Much of the literature on this topic focusses on 
the case of transitions between two vacua having different values of the cosmological 
constant. In this special semiclassical calculation of the nucleation rate based 

on instanton methods has been presented by Coleman and De Luccia [6]. However, this 
calculation, as well as many alternative approaches developed by others (cf. e.g. [37]), 
heavily relies on a high degree of symmetry of spacetime, which is initially assumed to be 
in a pure vacuum state with the geometry of de Sitter spacetime. We therefore feel that 
the applicability of these results to cases where spacetime is not (or only approximately) 
in a vacuum state is in need of some clarification (see also the discussion in [38]). 

A calculation of nucleation rates in arbitrary non-vacuum states, including all 
possible effects, is clearly beyond our capabilities. We will instead only take a small step 
away from the assumption of de Sitter symmetry and consider Friedmann universes in 
general, of which de Sitter spacetime is only a special case. We therefore retain many 
useful simplifications, in particular the assumption of homogeneity of space (but not of 
spacetime!), and the so-called 'thin- wall approximation'. The cosmological expansion 
of the universe, however, follows a non-trivial dynamical law (the Friedmann equation), 
and we are interested in the effect of the time-dependent Hubble rate on the nucleation 
rate of bubbles, which will in turn itself become time- dependent. 

To keep matters as simple as possible, we will consider the nucleation of a spherical 
bubble of new phase, and we will assume that its shell - the layer which separates 
the new phase from the old - is of negligible thickness compared to the size of the 
bubble. This amounts to the 'thin- wall approximation'. The energy budget of the bubble 
consists of latent heat (the difference between the energy densities of the two phases) 
and surface tension. Throughout this section we will neglect gravitational backreaction 
of the bubble onto the spacetime geometry since we are primarily interested in the 
effects of cosmological expansion of the background, which is taken into account, and 
the treatment of the full gravitational problem would introduce too many additional 
complications. 

2.1. Lagrangean formulation 

Associated to the background spacetime is the Friedmann-Lemaitre-Robertson- Walker 
(FLRW) hue element 

ds^ = a\r]) [-dr/2 + dr^ + r^Q^] , (2) 

where we have assumed flat spatial sections. The conformal time rj is related to 
cosmological (proper) time t by dt = adt]. 

Spherical symmetry reduces the bubble dynamics to a 1 + 1 dimensional problem. 
Denoting the coordinate radius of the shell as f, the shell trajectory f (77) follows from 
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the action 



S 



dr] 



'Air 



(3) 



where e and a denote, respectively, the difference between the energy densities of the 
two phases (latent heat) and the surface energy density (surface tension) of the shell. As 
we have indicated above, this effective action does not take into account gravitational 
self-interaction of the bubble. In order to include some of these effects, one could add 
surface-surface, volume- volume, as well as surface-volume terms for gravitational energy. 

The evolution of the scale factor a introduces an explicit time-dependence, giving 
rise to nucleation rates which will in general be time-dependent as well. A formalism for 
calculating semiclassical tunneling rates in time-dependent settings has been presented 
by Keski-Vakkuri and Kraus [39]. Its application to the present scenario will be worked 
out in detail in the following section. 



2.2. The complex time path formalism 

Our starting point is the classical equation of motion of the bubble, which can be found 
from eq. (3) as 



47reaV^ - STraaVWl - (d^fY = — 
^ drj 



1 - {dnV 



-\2 



(4) 



and we have dropped some labels in favor of simplified notation. 

The classical trajectory after tunneling emanates from a classical turning point 
where the canonical momentum 
dC 



P 



ddrfT 



Ana- 



(5) 



vanishes. It has been pointed out in [39] that by analytic continuation to complex t] 
one can find a classical trajectory|| (in the complex rj plane) that smoothly shrinks the 
bubble to zero size. To this end, it is useful to rewrite eq. (4) as an equation for rj (f): 



Airea r dfrj — Siraa ry {dfrj) 



dr 



Ana 



which, after some simplification, becomes 

lav/(5^ = 2^ + 3^ 
a ^ r a 



1 



(6) 



(7) 



II It turns out to be a special feature of time-dependent settings that this trajectory is not along a purely 
imaginary direction as would be the case in static settings, where, as a consequence, Euclideanization 
of time is a valid prescription. 
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We are looking for the solution with the boundary conditions 



p{r]o) = Ana 



, drV (0) = . (8) 



ri=VO 

The first condition matches the solution to the turning point, from which on it remains 
on the real rj axis for increasing f. The second condition guarantees that the full 
spherically symmetric solution is regular at the origin f = 0. Since r) will depart from 
the real axis for radii smaller than the nucleation radius, it is clear that one also has 
to analytically continue the time-dependent scale factor a to the complex plane. These 
boundary conditions do not determine the solution entirely: we are still free to choose 
the 'nucleation time' rjo. The 'nucleation radius', i.e. the coordinate radius of the bubble 
at the classical turning point, will accordingly be denoted as fo, and fulfills ?7 (fo) =770. 
This means that we have a one-parameter family of solutions labeled by their individual 
nucleation times. For each solution, the semiclassical tunneling rate is determined by 
the imaginary part of its action: 

r(r/o) ~exp[-2Im5(r/o)] , (9) 
where it is again useful to write ImiS (?7o) in terms of 77 (r/o! ^)'- 



ImS (?7o) = Im / df 
Jo 



^ea'^{r] (770; r)) r^dfi] (770; r) 



-Airaa^ir] (770; r)) r^y^ [dfT] {r]o;r)f - 1 



(10) 



We have now all the tools to compute time- dependent tunneling rates. Let us now 
turn to some explicit examples. 

2.3. Minkowski spacetime 

Before dealing with time-dependent backgrounds, let us briefly review the situation 
in Minkowski spacetime. We can simply set a = 1 and rj = t. Since we know that 
the solution has to be invariant under boosts, a natural ansatz for the trajectory is a 
hyperbola 

-{v- = rl , (11) 

which yields 

viVo^r) =7/0 + -fg , (12) 

where the positive sign is chosen for the square root, corresponding to an expanding 
bubble for f > fo. Inserting into eq. (7) one infers 

ro = f , (13) 

and one can check that this solution fulfills the boundary conditions (8). The imaginary 
part of its action can be found readily from eq. (10) , 

Im^Mink = —erg = , (14) 
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which is exactly Coleman's result [9]. 
2.4- de Sitter spacetime 

Our first example of an expanding universe will be de Sitter spacetime. Although it can 
actually be written in static coordinates, the FLRW metric with flat spatial sections has 
a scale factor which grows exponentially with time, a = exp (Ht). Written in conformal 
time this becomes a = —l/Ht], where rj runs from — oo to as t runs from — oo to +oo. 
In this case, eq. (7) reads 



\2 -, c^dff] 3 djr] 



„ {drvY -1 = 2^-- • (15) 

While finding the complete solution to this problem does not seem promising, it 
turns out that we can guess the relevant solution by sensibly generalizing eq. (12). Since 
de Sitter spacetime (in flat coordinates) has a constant expansion rate H, we expect 
that the proper nucleation radius should be independent of time. Now, if fo is to be 
a comoving radius, it has to be divided by the scale factor, i.e. instead of eq. (13) we 
expect 

_ _i / X 3cr 3cr 

ro = a (?7o) — = -HrjQ— . (16) 

Remarkably, choosing this nucleation radius is enough to have eq. (12) solve eq. (15) 
with the boundary conditions (8). 

The integral of eq. (10) is still solvable, and its imaginary part is found to be 

2 



Im5ds 



I - \ll + {?>Ha/ef 



sjl + {?,Ha/ef 

= ^sinh2^1n(l + (3i7a/6)^) , (17) 

which is independent of the choice ?7o for the time of nucleation. This means that the 
nucleation rate in de Sitter spacetime is time independent, which is a manifestation of 
the fact that de Sitter spacetime has no true dynamics. This expression reduces to the 
result of Minkowski spacetime eq. (14) in the limit oi H ^ 0. The first correction is of 
order and complies with the expansion obtained by Abbott, Harari and Park [40]. 
We can also take the limit e ^ 0, which corresponds to the nucleation of a domain wall 
separating two degenerate vacua. One finds 

limlm^ds = , (18) 
in complete agreement with a result obtained by Basu, Guth and Vilenkin [41]. 

2.5. More general FLRW spacetimes 



Thus far we have only considered static spacetimes in order to get some experience using 
the new tools. Let us finally turn to more general spacetimes, where the expansion 
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rate is not assumed to be constant. A simple deformation of de Sitter expansion is 
given by power law inflation, where the scale factor grows as a = (771/77)^^°. For small 
deformation parameters a this is an exact slow-roll solution of inflation with a constant 
slow-roll parameter —dtH/H"^ ^ a. rji denotes an (arbitrary) point in time where the 
scale factor is normalized to unity. This spacetime is obviously not static, and we shall 
use this simple example to study effects of time-dependent cosmological expansion on 
tunneling rates. For power law inflation, eq. (7) reads 



We have been unable to find an analytic solution to this equation and therefore decided 
to treat it numerically. Using the boundary conditions (8) one can find numerical 
solutions for any choices of a, rjo/rji and e/a. A parameter study of the tunneling rates 
reveals the following picture. There are now three different time scales in the problem. 
The inverse expansion rate H^^ (rj) gives the time scale on which the background (scale 
factor) changes significantly. But there is now another time scale related to the change 
of the expansion rate itself (higher order derivatives of the expansion rate are zero in 
power law inflation). These two time scales have to be compared to the 'bubble crossing 
time', which we define by the nucleation radius divided by the speed of light, or roughly 
3a /e in our units. 

If the bubble crossing time is the smallest time scale of the problem, the tunneling 
rate is well approximated by the result of Minkowski spacetime, eq. (14). However, if the 
bubble crossing time is not much smaller than the Hubble time at nucleation, if ~^ (770), 
there are two different possibilities. Either the characteristic time scale for the change of 
the expansion rate, \dtH/ H\~^ , is still much larger than the bubble crossing time - then 
we are in a regime where a quasistatic approximation is valid such that a good estimate 
for the tunneling rate can be obtained from eq. (17) by setting H = H^rjo). Or the 
bubble crossing time cannot be regarded as small with respect to any other time scale. 
In this case, the tunneling process 'feels' the changing expansion rate, and the decay 
rate is modified significantly. Our numerical study clearly indicates that the tunneling 
rate is enhanced compared to the quasistatic estimate. We believe that this is related to 
the fact that the expansion rate decreases with time. Instead of setting H = H (rjo) in 
eq. (17), one should average the expansion rate over an interval of one bubble crossing 
time prior to nucleation. Using this averaged expansion rate in eq. (17) turns out to 
give a much better estimate for the actual tunneling rate, cf. Fig. 2. 

2.6. Spacetimes with Big Bang singularity 

An interesting case which is also relevant for the chain inflation scenario is a universe 
filled with radiation and some vacuum energy, such that radiation dominates its early 
evolution, while it approaches a vacuum de Sitter phase at later time. If we neglect the 
brief era of matter domination, this is also a good model for our present universe (after 
reheating has taken place). During the radiation era, the Hubble rate drops rapidly. 




(19) 
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Figure 2. Numerical values for the 
imaginary part of the action S for 
tunneling solutions with proper nucleation 
radius Scr/e = 2H^^{riQ), as a function 
of the slow roll parameter dtH/H^. The 
values are normalized to the 'quasistatic 
approximation', which is obtained from 
eq. (17) by setting H = H{r]o). Thus, 
in this approximation only instantaneous 
dynamical parameters are taken into 
account. The tunneling process, however, 
has a characteristic time scale given by the 
light-travel time across the bubble. Hence 
it is more appropriate to take into account 
some average dynamics of the background. 
In a crude way, this can be accomplished 
by using an averaged expansion rate to 
evaluate the tunneling probability in the 
quasistatic approximation. The dashed 
red line is obtained from a proper time 
average of H taken over an interval At = 
3a /e before nucleation, whereas the dotted 
green line is obtained from a conformal 
time average taken over an interval Ar/ = 
At/a{r]o). 

and assuming the vacuum is metastable, one might wonder if this affects tunnehng 
rates as described above. We can answer this question by comparing the relevant time 
scales. If the universe is sufficiently flat, a oc t^^^ as long as radiation dominates over 
vacuum energy. Note that the universe has a Big Bang singularity at t = as long as 
reheating and any kind of cosmology before the radiation era is not taken into account. 
Hence, there exists a particle horizon of size ~ t. Also both the Hubble time and the 
characteristic time scale for the change of the Hubble rate are of order ~ t. Since even 
tunneling cannot violate causality, bubbles larger than the particle horizon cannot be 
produced, at least not in the semiclassical picture we are using here. This means that 
the bubble crossing time can never be much larger than the other relevant time scales, 
and thus we expect corrections to the tunneling rates to be small in general. Moreover, 
a numerical study indicates that the tunneling rates for bubbles whose nucleation radius 
is comparable to the size of the particle horizon become sensitive to the dynamics of the 
scale factor up to the vicinity of the Big Bang, such that details of reheating and even 
earlier cosmology become relevant. 

As soon as vacuum energy begins to dominate, the time scale on which the 
expansion rate changes tends to infinity quite rapidly as the Hubble rate approaches a 
constant. Hence, for the late part of the evolution, we expect the de Sitter approximation 
to be good. Note, however, that a particle horizon still exists and that bubbles should 
not exceed the horizon size at any time as long as causal physics is at play. To eliminate 
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this 'horizon problem'^, one has to change the details of the model near the Big Bang 
singularity to include any cosmology preceding the era of radiation domination. 

3. Bubble propagation on dynamical backgrounds 

In this section we will explore the effect of background irregularities on the evolution 
of the bubble and whether these effects can potentially be seen by an observer inside 
the bubble. To this end we will follow a somewhat different route than in the preceding 
section, where we ignored gravitational backreaction of the bubble. In order to include 
this backreaction, we will follow the approach used in [29] and assume that the bubble 
wall separates spacetime into two parts, described by different metrics and containing 
different matter. These two parts, for convenience called interior and exterior part, are 
approximated by the manifolds A4- and Ai+ and are joined along a common timelike, 
spherically symmetric hypersurface S which represents the bubble wall. Since the work 
of Israel [42] there exists a well established formalism for joining manifolds along a 
common boundary, known as the Israel junction conditions. Once Ai- and Ai^ are 
given, the evolution of the bubble wall uniquely follows from these conditions. Moreover, 
by this construction, the resulting spacetime is a solution to Einstein's field equations. 

For the interior part we use de Sitter spacetime (which can be thought of as an 
approximation to the inflationary phase of our observable patch of the universe) and 
for the exterior part two different cases will be considered. We will study the evolution 
of vacuum bubbles in an inhomogeneous background for which the exterior part is 
approximated by the spherically symmetric Lemaitre-Tolman-Bondi (LTB) spacetime 
[43, 44, 45]. To maintain spherical symmetry it is assumed that the bubble nucleates 
in the center of the LTB model. In a similar vein, we will also explore the evolution 
of vacuum bubbles when the exterior part is given by a FLRW universe filled with a 
fluid which at some time undergoes a smooth (second-order) phase transition, e.g. from 
w = —1 to w = 1/3, or vice versa. 

In the first subsection we will introduce the description of the bubble wall and the 
interior and exterior parts of spacetime. Thereafter we will give a concise guide on how 
to calculate the junction equations and write them down explicitly for the cases of our 
interest. Finally we will solve these equations numerically and discuss the results. 

3.1. Bubble wall and background spacetime 

3.1.1. Interior: de Sitter The interior of the bubble is assumed to be in a de Sitter 
phase with vacuum energy density given by A/ (Svr). We employ the flat slicing in which 
the metric is given by 



^ Actually, this is nothing but the good old horizon problem of standard Big Bang cosmology in a new 
guise, and therefore it can be solved in the same spirit. 




(20) 
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and stress-energy is given by T^^, = —Ag^^. All quantities should carry the index (— ) 
which indicates that they belong to A^_. For convenience we write these indices only 
when necessary, hoping that it will be clear from the context which quantities are meant. 

3.1.2. Bubble wall The timelike, spherically symmetric hypersurface separating the 
interior and exterior parts of spacetime, the bubble wall, is described by the metric 

ds^ = -dr^ + i?2dfi2 . (21) 

Though it has been shown [46] that the stress-energy of a wall which separates regions 
of different vacua will solely be given by the surface tension, it is not clear whether 
this conclusion is valid when matter is present. However, in lack of a field theoretic 
description, we assume that stress-energy is given by 

Sij = —crhij . (22) 

where hij is the metric tensor of (21). The equations of motion of the proper radius R 
and surface tension a of the bubble will be given by the junction conditions. 



3.1.3. Exterior: LTB and FLRW 



LTB spacetime In order to explore the evolution of a bubble in an inhomogeneous 
background we use an LTB ansatz. In comoving coordinates a suitable metric [47] can 
be given in the form 

d.- = -if + '-"' dr' + a\t, r)rW , (23) 

1 -|- 2h[r) 

with 2E{r) > —1 but otherwise arbitrary. From Einstein's equations with a dust source 
Tf_tu = pSji^i — ^Qfiu , we obtain the equation of motion of the scale factor 

Here, M(r) is the first integral of motion which corresponds to the active gravitational 
mass within a sphere of coordinate radius r. Once the scale factor is known, the dust 
density p is determined by 

Svrp - ^^'■^ (25) 

g2j.2 (^^Q^d _|_ a) 

FLRW spacetime In addition, we want to study the motion of bubbles in a background 
which undergoes a phase transition. Therefore a flat FLRW spacetime is employed, a 
comoving coordinate system of which is 

ds^ = -dt^ + a^it) (dr^ + rMfi^) , (26) 
with stress-energy given by 

T^,. = {p + p) Sl6l + pg^, - Ag^, . (27) 
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The evolution of this background follows from the Friedmann equation 



^V = -P+-, (28) 
a y 3 3 ' ^ ' 

and continuity equation 

d,p + 3^{p + p) = 0. (29) 
a 

To look at the influence of a phase transition in the background on the evolution of the 
bubble, we artificially"'" introduce an abrupt change in the equation of state p = wp via 

w{t) = (1 ± 2tanh(7pt(t - V))) . (30) 

to model a nearly instantaneous (on time scale 7"^^ <^ H^^) phase transition at t = tpt 
from w = —1 to w = 1/3 ('reheating') and vice versa. Solutions to equations (28) 
and (29) can not be given in closed form and will be obtained numerically. 

3.2. Conditions for a valid junction 

The problem of joining manifolds across a common boundary has been lucidly explained 
in the work of Israel [42]. For a nice introduction see also the textbooks [48, 49]. Two 
conditions arise in the course of joining two manifolds. The first junction condition 
requires the induced metrics hij = g^u^^e'^ to coincide on S 

[hi,] = h^\,:-hr^\j, = . (31) 

The second junction states that, whenever there is a discontinuity in the extrinsic 
curvature of E as seen from a. surface layer of stress-energy Sij, given by 

SnS^j = [Kij] - hij [K] (32) 

will be present. Therefore the proposed stress-energy on the bubble wall (22) has to be 
identified with the difference in the extrinsic curvature. The components of the extrinsic 
curvature tensor Kij are defined as the covariant derivative of the vector e^ along e^ 
projected onto the surface normal 

= n^r^^efe,^ . (33) 

Once the projectors ef = dx^/dy^ are known, the normal vector of S can be obtained 
by the conditions 

n^n^" = 1 and n^ef = (34) 

up to a sign which determines how and are stuck together. We choose this 
sign such that 

n^, = y/g^ {-r,i, 0,0) , (35) 

^ A universe with two or more components with different equations of state, Hke e.g. radiation and 
cosmological constant, actually has one or several intrinsic phase transitions. However, these transitions 
are very gentle and would probably only produce a minuscule effect. 
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where a dot refers to a partial derivative with respect to r. This choice imphes that 
in radii increase towards S and in Ai^, radii decrease towards S. The continuity 
equation 







(36) 



is not independent of the two equations resuhing from (32) and will be used to substitute 
one of these. We consider exterior stress-energy given by a perfect fluid and interior by 
a cosmological constant. Therefore the r component provides the following first order 
equation for a 

a = (p + p) ^/g^ri . (37) 
In the next subsections we will write down these equations for LTB spacetime. 



3.2.1. Equation of motion for the size and surface tension of the bubble By virtue 
of spherical symmetry of the bubble wall, angular coordinates can be identified and 
the time and radial coordinate can be parameterized by the proper time of the bubble 
(t(r),r(r)). We write down the equations for the LTB part and relegate the FLRW 
equations to the appendix. The conditions of the first junction turn out to be 



ar 



R, = 1 + 



(3^ 



1 + 2E 

We will make use of the 6'^-component of the second junction condition which yields 



AnaR 



R^ 



A. 



i?2 + l 



2M A, 



3 "~ V ' ' " R 3 
Solving for the derivative we obtain the more convenient form 



-R' + V 



1 
2 



with 2V given by 



2V 



A. 



A+ - A. 
247ra 



2TTa 



R' 



1 



A+ - A_ 
487rV2 



M 

'r 



M2 



167rV2i?4 



(39) 



(40) 



(41) 



If M = const all coefficients in the potential are constant and we recover the 
Schwarzschild-de Sitter model which has been discussed in [46, 50, 51, 52]. 

However, since we want to introduce an exterior matter density, M will no longer be 
constant, i.e. M(r) = M{R/a) and the scale factor of the ambient spacetime will enter 
the equation. In this way the motion of the surface becomes sensitive to the presence 
of matter in the background. 

Note that (23) is covariant under a rescaling of the radial coordinate. Together 
with drM > this allows one to define a radial coordinate such that M(r) = ^Ar^ 
where A is a constant. The potential becomes 

2" 

2na 1 R^ . (42) 



2V 



A_ 
~3~ 



A 



A, 



A_ 



3a-V 



2A7Ta 
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and the equation of motion for the surface tension is 

& = p / rt . 43 
VI + 2E 

It is restricted by equation (39) to 



4.a<^^ + ^^^^-e. (44) 

Here e is the difference in energy density between inside and outside. This bound is a 
direct consequence of the geometry that was fixed by the sign of the normal vector (35). 

3.2.2. Evolution equations expressed in exterior coordinates Since the background 
dynamics of the LTB (and FLRW) part of the spacetime can be obtained only 
numerically all equations will be solved in these coordinates in the first place. Making 
use of the junction conditions we are able to write down the evolution equations in terms 
of the exterior coordinates. After a solution to these equations has been obtained, the 
matching conditions will be employed again to express the evolution in terms of the 
interior coordinates. Like before, we explicitly write down the expressions for the LTB 
part and provide the FLRW equations in the appendix. 

Let f{t) be the bubble radius in these coordinates. Writing R = i-^ (af) we can 
solve for dtf and obtain 



-(1 + 2E)fdta + X {1 + 2E) (1 + 2V) Urdtaf -2E + 2V) 

dff = ^- . (45) 

* {fdfa + a) {2E - 2V) ^ ' 

where we have chosen a positive sign of the square root because we are interested in 
solutions of physically growing bubbles, i.e. i? > 0. The equation for a (37) can as well 
be converted to LTB coordinates 

{fdra + a) dtf 

Ota = p I . (46) 

^Jl + 2E-{fdra + af{dtff 

The surface tension becomes time dependent. It increases when the comoving radius of 
the bubble increases, i.e. it collects matter from the background, and if the bubble 
shrinks it will exactly provide the amount of matter density determined by the 
background. This is a limitation of the spacetime junction approach, which in the 
present form does not capture the physics of matter transfer across the junction surface. 
Physically we would expect that dust would actually penetrate into the bubble, as one 
can convince oneself by looking at the geodesies of 'test particles'. However, interior 
and exterior parts of spacetime are fixed ab initio and cannot be changed by the motion 
of the bubble. Bearing with this limitation, we continue our analysis and will give an 
outlook on this issue in our conclusions. 

For now, equations (45) and (46) completely determine the evolution of the bubble 
in exterior coordinates. 
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3.3. Bubble evolution on dynamical backgrounds 

In this section we explore the propagation of bubbles on dynamical backgrounds. We 
start with an exact solution of a de Sitter/de Sitter spacetime and continue with the 
numerical solutions obtained for the de Sitter/LTB and de Sitter/FLRW spacetimes. 



3.3.1. Exact solution in de Sitter/de Sitter spacetime We solve equations (45) and (46) 
for the case that both spacetimes are de Sitter with cosmological constants given by K±. 
Then M = and we employ the flat slicing where also E = {]. Since there is no matter in 
the background it follows from equation (46) that a = const. The potential V reduces 
to 



2V 



For further convenience we define 



. (47) 



(4^ 



where H'^ 



A-I-/3. The solution is valid in 
[u_,A_) and (n_|.,A+) respectively. 



such that 2V = -Hl{l + u\)R^ 
M.^ and A^+, with the corresponding quantities 
Equation (45) becomes 

?}fM^.l_HV_2V.(,,)_l, (49) 
H r 

which, when rewritten as a differential equation for V , can be solved by separation of 
variables. Solving for f yields 



fit) = ^^-2 + (exp i-Ht) - ifH-^ . 

For convenience we normalized the scale factor at the time t 
took to = without loss of generality. This implies 

e 



to when dtr{tQ 



± 



3a 



=F 27ra 



(50) 
and 

(51) 



After t = to the bubble accelerates and converges to r{t oo) = yl + u^ro, see Fig. (3). 
We conclude with the remark that (50) reduces to the solution found in (16) in the limit 
of G ^ 0. 



3.3.2. Numerical solution in de Sitter/LTB spacetime The goal of this section is to 
understand the influence of ambient inhomogeneities on the motion of the bubble wall. 
The first step into the numerics of the bubble is solving the dynamics of the background. 

For an intuitive approach we define 2E{r) = —k{r)r^ where the profile k may 
be interpreted as the local spatial curvature. Then, the coordinate size of the spatial 
section, if finite, is determined by k (r^ax) '^max = 1- In addition, we have to specify the 
initial value of the scale factor ao(^) = a(to, r) which will in general depend on the radial 
coordinate. Instead, one may equivalently choose the initial dust density po{r) = p(to, r), 
which defines ao{r) via equation (25). Thus, in coordinates where M(r) = ^Ar^, spatial 
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Figure 3. Trajectories of the bubble wall 
in the flat slicing of de Sitter/de Sitter 
spacetime for different values of u ; the pa- 
rameter u (48) encodes the dependency on 
A± and the surface tension a. The bub- 
ble expands accelerated from coordinate 
radius tq — u^^H"^ but converges to the 



finite coordinate radius vl-M?fo in the 
limit i — > oo. The trajectory is shown in 
de Sitter coordinates rather than in physi- 
cal quantities (proper radius R vs. proper 
time r) to make comparison to our later 
results easier. 
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inhomogeneity of the LTB spacetime is incorporated in the functions k and po- In the 
Umit where k and po a^rc constant the model becomes homogeneous. 

After both functions have been specified, the partial differential equation (24) can 
directly be integrated at each r. When a solution is found its validity has to be checked. 
If not dr{ar) > 0, the weak energy condition p > is violated by the occurrence of 
a shell-crossing singularity. This actually restricts the curvature profile k, since, if too 
steep, it disturbs the background massively such that the dust density will violate the 
weak energy condition at some time. We evolve the system until the background space 
has expanded for about 4 efolds. 

The initial time of the analysis is t = to where dtf{to) = 0. to will be referred to 
as the time of nucleation of the bubble. The nucleation radius is determined by the 
parameters A, A+, A_, ao. It can be inferred from equation (45) which we rearrange to 



Note that the initial difference in energy density eo includes the initial dust density po, 
which can be a function of fo, too. This equation illustrates the route that we will follow. 
There are two possibilities, via the functions k and po, to introduce inhomogeneity in 
the LTB model. These two cases are considered independently, meaning that one of the 
two terms on the right hand side will be independent of r. In case there is more than 
one solution to the equation the smallest positive value will be taken. 

Note also that the proper kinetic energy of the bubble at nucleation is proportional 
to if\t=to = H^Rl, where Ho = ^{to,fo) and Rq = ao(fo)fo. 

Homogeneous limit The first thing we will explore is not the effect of inhomogeneity, 
but simply what happens when the bubble nucleates in a background where, in addition 




(52) 
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to vacuum energy density, some dust is present. To keep things simple we refrain from 
adding any curvature at this stage and set k = and ao(r) = 1. The spatial sections of 
LTB spacetime become homogeneous and reduce to the FLRW limit. Henceforth, we 
will always assume that at nucleation dust density shall exceed exterior vacuum energy, 
and for definiteness we choose SttA/S = 10~^ and A+/3 = 10~^ in accordance with [29]. 

It is important to note that this setup already has a significant effect on the 
evolution of the bubble. Whereas a bubble that nucleates in vacuum always begins 
to expand, this is no longer guaranteed as soon as considerable amounts of matter are 
present in the environment. This can be seen by the following argument. The force 
which accelerates the shell has two contributions: one from the surface tension, which 
is always directed inwards, and one from the pressure difference between the interior 
and exterior fluid. In the case where both fluids are mere cosmological constants it is 
easy to show that the pressure force can sustain the surface tension and will push the 
shell outwards. However, if the exterior fluid is mainly composed of pressureless dust, 
at some point surface tension will outrun pressure support and the bubble will be forced 
to collapse. 

To make this statement more quantitative we can take another derivative of eq. (45) 
and examine the behavior of the bubble when dtf = 0. In the homogeneous limit one 
finds 

t \dtr=o a \^ 247ra 3a J ^ ^ 

The sign of dff depends on how o"^ and p compare to the latent heat of the vacuum, 
evac = (A+ — A_) /Stt. The bubble can only expand into the ambient spacetime if 
p is not too large. In particular, one can never have an expanding bubble during a 
matter dominated phase* . This strongly limits the possibility to study the propagation 
of bubbles into inhomogeneous matter with the current approach, since the exterior 
spacetime has to be vacuum dominated in order to allow the bubble to propagate 
towards the inhomogeneity in the first place. We have summarized the behavior of 
freshly nucleated bubbles of vacuum within a dust environment in Fig. 4. 

We numerically calculated the trajectories of bubbles in a dust dominated 
background. As expected by the argument above, in contrast to the de Sitter case, 
bubbles contract as seen from the exterior perspective (Fig. 5, left). In fact some 
bubbles contract so fast that the growth of their physical size is decelerated, i.e. their 
proper kinetic energy decreases. Our results show that it will even decrease to zero for 
small bubbles with HqTo < 1. In this case the proper kinetic energy becomes imaginary 
and the evolution of the bubble had to be stopped. When fo > Hq^ they retain some 
proper kinetic energy but nevertheless shrink in exterior coordinates and converge to a 
coordinate radius which is smaller than the coordinate radius of nucleation. For bubbles 

* Note that this was not at aU an issue in section 2, since background spacetime was assumed spatially 
homogeneous and any dust would therefore permeate the bubble. In the present setup, however, the 
interior is assumed to be completely empty except for a possible cosmological constant. In this sense, 
the issue is one of initial conditions. 
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Figure 4. This plot characterizes the 
early-time behavior of a vacuum bubble 
after it was nucleated at rest in the 
comoving frame of an exterior flat FLRW 
spacetime with dust and cosmological 
constant. With respect to an exterior 
comoving observer, the bubble shows 
different behavior in different regions of 
the cr^-/9-planc, which is drawn in units of 
Evac = (A+ — A_) /Stt. The shaded region 
is forbidden for our choice of junction 
because a there violates the bound (44). 
If the energy density of dust p is chosen 
above the dashed red line, the bubble 
starts to contract. This includes all matter 
dominated universes p > pvac = ^+ /Stt 
since we assume A+ > A_ > 0, which 
implies pvac > Cvac- Below the dashed 
red line the bubble starts to expand into 
the ambient spacetime. This includes all 
vacuum dc Sitter spacetimes since they 
are found on the line p = 0. Below the 
dotted blue line, all universes are vacuum 
dominated. 



larger than 2Hq^, even if they fulfill the bound (44) initially, the dust density in the 
background drops faster than the surface tension of the bubble such that the bound will 
be violated soon. 

Note also that after a/A_|./3 icq — 0.4 the bubble moves on a vacuum dominated 
background. The contour plot in 5 illustrates the fate of a bubble in dependence of the 
parameters (Aq, o"o). 

Inhomogeneous dust density We now explore possible effects of inhomogeneity by the 
introduction of an initial dust distribution po{r). The bubble may nucleate either in an 
over dense or in an underdense region of space. It is not very revealing to consider a 
radially decreasing dust profile, because the dust will be rarefied by the expansion of the 
background anyway. Therefore we will have a look at radially increasing dust profiles 
only. For a given pQ{r) , the initial scale factor is determined by 



K.2 



al{r) = ^41 ^dr. (54) 
H J po{r) 

We will consider the profile 

Po(r) = Aryr\ (55) 

— 1 /2 

with rA = (A4./3) ' . We expect that small bubbles with fo -C will begin to expand 
initially because they will find themselves in a vacuum dominated background, whereas 
large bubbles with fo > ta are in a matter dominated background with their subsequent 
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Figure 5. Results obtained for tlie liomogeneous limit of the LTB model with k ~ 
and ao(r) = 1. Left figure: Trajectories of the bubble wall for several values of the 
surface tension Q.SHq < Atto-q < 0.75i?o and A_ = 0.1A+. The fate of the bubble 
depends on the nucleation size. Small bubbles contract in LTB coordinates until their 
proper kinetic energy becomes zero, i.e. R = and the equation of motion becomes 
imaginary. The greater the bubble the more likely it sustains kinetic energy until the 
background is dominated by A_|_ and it converges to some finite coordinate radius. 
Right figure: Fate of a bubble in dependence of the parameters ctq, A_. We looked 
in the region of parameter space where the nucleation radius of the bubble is within 
0.1 < Hofo < 5. The black lines are lines of constant kinetic energy and nucleation 
radius Hofo = (0.1, 1, 2, 3, 4, 5). The black shaded region is considered to be unphysical 
since equation (44) is not fuUfiUed there already initially. In the region shaded light 
blue, bubbles either contract until 7? = 0, or they hit the geometrical bound (44). 
Neither occurs in the white region in which bubbles 'survive' and come to rest at a 
finite coordinate radius. 



evolution being much like in the homogeneous limit discussed before. Although the 
initial dust profile is increasing, expanding bubbles propagate into regions of lower 
density due to the expansion of the background, see Fig. 6. 

Inhomogeneous curvature The other possibility is to incorporate inhomogeneity in 
the neighborhood of the bubble via a curvature profile k. However, in view of the 
results obtained in the homogeneous limit, we state that the bubble will hardly be able 
propagate into that inhomogeneity because it will shrink as long as the background 
is matter dominated. Even if the condition that the bubble nucleates comovingly is 
relaxed, such that the bubble may have dtf(to) > 0, deceleration is large enough to 
make the bubble contract almost immediately. Therefore it seems that studying the 
effect of curvature inhomogeneity within this approach is hardly feasible. Note that this 
result appears to be in contrast to what has been obtained in [29]. 

Nevertheless, to see what happens when a bubble enters a curvature inhomogeneity 
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Figure 6. Results obtained for the inhomogencous LTB model with the dust 
profile (55). The left figure shows the trajectory of bubbles with initial surface tension 
0.7\/A+/3 < 47rcro < 1.6a/A+/3 and A_ = 0.1A+. Smaller bubbles nucleate in 
a region where vacuum energy dominates over dust density and therefore expand. 
Bubbles with fo ~ contract because they are already in a dust dominated 
background. The right figure shows the dust density at the position of the bubble. 
Although the dust profile (55) radially increases, the expansion of the background 
dilutes matter efficiently such that expanding bubbles effectively propagate in a 
decreasing profile. 



we make use of a result obtained previously. In the last section it was shown that bubbles 
which nucleated in a vacuum dominated region expanded initially. Now, we add some 
curvature in 'front' of the bubble and have a look what happens if the bubble encounters 
that inhomogeneity and whether there is a difference compared to the corresponding 
solution in the flat background. The curvature profile is given by 

where TZct = {4:TrAy/A^)~^^^ . Taking the dust profile from the last section, we fix the 
initial surface tension and vacuum energy of the bubble to (To = 0.6^A+/3, A_ = 0.1A+. 
The evolution of the bubble is significantly affected in the exterior perspective. However, 
this may be just a coordinate effect and when looked at the trajectories in the interior 
coordinates the effect practically vanishes. Nevertheless, there remains an effect in the 
surface tension of the bubble, see Fig. 7. 

3.3.3. Numerical solution in de Sitter/FLRW spacetime In this section, the evolution 
of the bubble on a de Sitter/FLRW background will be considered. The FLRW part 
is supposed to contain vacuum energy and a perfect fluid which undergoes a phase 
transition i) from w = 1/3 to w = —1, or ii) vice versa. The FLRW dynamics (28) 
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Figure 7. Results obtained for a bubble propagating in an inhomogcneous background 
(curvature inhomogencity and inhomogeneity in the initial dust profile). Upper left: 
the curvature profile (56) with ai = 1.5, 012 = 0.02 and ~ 0.2. Lower left: 
Bubble trajectories in the comoving LTB coordinates in a flat background (red) 
compared to the curved background (green). The trajectories are affected significantly 
in these coordinates. However, if one converts the trajectories to the interior de Sitter 
coordinates the effect practically vanishes (upper right). Nevertheless, there remains 
an effect on the surface tension of the bubble (lower right). 



and (29) will be solved numerically with the initial values a(to) = 1, 87rpo/3 = 10"'' 
and A_|_/3 = 10~^. Again = shall represent the exterior time coordinate at 
which dtf{to) = 0. The phase transition is supposed to occur at tpt = 0.5Hq^, and 
we set the width 7"^^ = 1 <^ Hq^ in order to model a nearly instantaneous transition. 
After these dynamics have been established we consider the evolution of the bubble. 

The nucleation radius of the bubble can still be obtained from equation (52) with 
k = 0. Of course, an immediate effect of the phase transition on the trajectory of the 
bubble can be seen in the exterior coordinates. In case i) the contraction of the bubble 
is stopped when the background becomes vacuum dominated, whereas in case ii) the 
expansion of the bubble reverses to contraction due to the phase transition. However, 
contrary to the inhomogeneous background discussed before, the effect is still present 
when the trajectory is expressed in the coordinates of an interior observer, see Fig. 8. 
However, the most prominent effect still can be seen in the surface tension of the bubble. 
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Figure 8. Evolution of bubbles in the de Sitter/FLRW spacctime for surface tensions 
OAHq < Ana < 0.6iJo and A_ = O.IA4.. The left column shows the results of a 
transition from w = 1/3 to w = —1 and the right column from w = — 1 to w; = 1/3 
compared, respectively, to their counterparts where the equation of state remains 
constant. The upper plots show the trajectory of the bubble in exterior coordinates, 
while the mid plots show the trajectories as seen from an interior observer. Unlike for 
the inhomogeneous background, we see that the exterior phase transition indeed leaves 
a sizeable effect on the trajectory of the bubble. In addition, the most prominent effect 
is again on the surface tension. 
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4. Conclusions 

First-order pliase transitions, as may liave occurred multiple times in the early universe, 
proceed by spontaneous nucleation of bubbles. In order to answer the question on how 
much information about the initial state 'survives' the phase transition, it is important 
to understand bubble nucleation and propagation beyond the assumption of a trivial 
initial state. This question may indeed be relevant for certain cosmological scenarios, 
such as the chain inflation proposal. In this particular scenario, a series of first-order 
phase transitions proceeds very quickly, such that the time in between two transitions is 
too short for the universe to dilute all inhomogeneities and thermal radiation produced 
in each transition. 

We have seen that the dynamics of the background spacetime affects the calculation 
of semiclassical tunneling probabilities. A simple comparison of time scales helps to 
decide if this effect is relevant. In cases where the bubble crossing time (the nucleation 
radius of the bubble divided by the speed of light) is much smaller than any time 
scale of the background geometry, we found that it is well justified to use tunneling 
rates obtained from field theory on Minkowski spacetime. Noting that this tunneling 
probability drops exponentially with increasing nucleation radius, it seems reasonable 
to assume that the bubbles in fact nucleate at tiny radii whenever we demand an 
appreciable tunneling rate which can lead to an onset of the phase transition before the 
universe reaches (approximately) a vacuum state. Hence, using Minkowski spacetime as 
an approximation appears self-consistent in many cases. However, if the expansion rate 
of the universe is very large, there may be cases where the effect of background evolution 
on tunneling rates can be important. In the context of the chain inflation proposal, we 
expect this to be the case very early on during a radiation dominated phase short after 
a previous phase transition. 

The existence of a particle horizon in a radiation dominated FLRW universe renders 
the nucleation of bubbles larger than this horizon impossible. This means that one has 
to have better knowledge of the spacetime near the singularity, including details about 
reheating and any cosmology preceding the radiation era, in order to avoid this 'horizon 
problem'. In other words, the pure radiation dominated FLRW universe is no useful 
approximation for a semiclassical calculation of nucleation rates of bubbles larger than 
the particle horizon, since those rates would be sensitive to the cosmology at the Big 
Bang. 

Concerning the subsequent evolution of comovingly nucleated vacuum bubbles in 
non-vacuum backgrounds, we have seen that already the presence of homogeneously 
distributed matter has a significant influence on the bubble. Unlike in de Sitter 
spacetime where the proper kinetic energy of a bubble increases exponentially, it may 
decrease in the presence of matter in the background. For small bubbles, the proper 
kinetic energy became zero and real classical solutions could not be obtained beyond 
this point. As has been pointed out in [29], such bubbles do not correspond to classical 
configurations, but should be interpreted as mere fluctuations. Bubbles with greater 
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radius have more proper kinetic energy and are able to survive until the matter density 
in the background has been diluted sufficiently. We note, however, that the setup in 
which we study bubble trajectories, in particular the choice of initial conditions, is 
somewhat ad hoc. To settle this issue, one would have to solve the tunneling problem 
with matter and gravity, an enterprise on which we did not embark in this work. 

We also studied the effect of inhomogeneities in the background on the propagation 
of vacuum bubbles. The results show that the trajectory of a bubble is affected, from 
the point of view of an exterior observer, when it propagates into an inhomogeneous 
background. Since it is not clear whether this is just a coordinate effect we have looked at 
the trajectory as seen by an interior observer. This observer will hardly see an influence 
of ambient inhomogeneities in the trajectory of the bubble. Nevertheless, when looked 
at from the inside, there remains an effect on the surface tension of the bubble. 

Furthermore, bubbles moving in an FLRW spacetime which itself undergoes a phase 
transition have been considered. Similar to the inhomogeneous case, a large effect is 
seen in the exterior coordinates. However, when converted to interior coordinates, there 
remains an appreciable effect in the trajectory of the bubble. This raises the question 
whether those perturbations of the bubble wall will lead to potentially observable effects. 
In the context of bubble collisions [53, 54], it has been pointed out that a disturbance 
in the trajectory of the bubble wall may lead to a redshift of the reheating surface and 
could therefore potentially be observable in the CMB. 

However, the most prominent of the effects of inhomogeneity or phase transitions in 
the background is found in the surface tension of the bubble. This is a consequence of the 
spacetime junction approach. As soon as a bubble propagates in a matter environment 
the evolution of its surface tension is determined by the demand of the background. 
This means that an expanding bubble 'collects' matter from the background while a 
contracting bubble will lose the amount of energy required by the space it uncovers. 
Therefore it is necessary to find a proper treatment of energy transfer through the 
bubble surface. Rather than fixing interior and exterior spacetime ab initio, one should 
dynamically construct those spacetimes from initial conditions. To solve this problem, 
an additional equation is needed, which arises from a proper dynamical description of 
the surface tension. It should capture the physics of matter transfer across the junction 
hypersurface and should probably be motivated from a field theoretical point of view. 
We hope to make progress in this direction in our future work. 

Acknowledgments 

We thank Ben Freivogel, I-Sheng Yang, Anthony Aguirre, and Douglas Spolyar for 
helpful discussions and suggestions. This work was supported by the German Research 
Foundation (DFG) through the Research Training Group 1147 'Theoretical Astrophysics 
and Particle Physics'. 



Tunneling and propagation of vacuum bubbles on dynamical backgrounds 



27 



Appendix A: Evolution equations in de Sitter/FLRW spacetime 

For the proposed FLRW spacetime the relevant equations can be obtained in the same 
way as for the LTB part by setting drtt = and E = 0. 

ar = R, = 1 + a^r^ . (A.l) 

The potential V reduces to 

2" 



2V = 

and the equations of motion are given by 
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, (A.2) 



fdta - J{1 + 2V) {{fdtay + 2V) 
d^f = ^ , (A.3) 

d,a=ip + p)-=£^=. (A.4) 
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These equations together with the background dynamics given by (28) and (29) 
determine the bubble motion in the de Sitter/FLRW background. 
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